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During the transition from a neutral to a fully reionized universe, scattering of cosmic microwave 
background (CMB) photons via free-electrons leads to a new anisotropy contribution to the temper- 
ature distribution. If the reionization process is inhomogeneous and patchy, the era of reionization 
is also visible via brightness temperature fluctuations in the redshifted 21 cm line emission from 
neutral Hydrogen. Since regions containing electrons and neutral Hydrogen are expected to trace 
the same underlying density field, the two are (anti) correlated and this is expected to be reflected 
in the anisotropy maps via a correlation between arcminute-scale CMB temperature and the 21 
cm background. In terms of the angular cross-power spectrum, unfortunately, this correlation is 
insignificant due to a geometric cancellation associated with second order CMB anisotropies. The 
same cross-correlation between ionized and neutral regions, however, can be studied using a bispec- 
trum involving large scale velocity field of ionized regions from the Doppler effect, arcminute scale 
CMB anisotropies during reionization, and the 21 cm background. While the geometric cancellation 
is partly avoided, the signal-to-noise ratio related to this bispectrum is reduced due to the large cos- 
mic variance related to velocity fluctuations traced by the Doppler effect. Unless the velocity field 
during reionization can be independently established, it is unlikely that the correlation information 
related to the relative distribution of ionized electrons and regions containing neutral Hydrogen can 
be obtained with a combined study involving CMB and 21 cm fluctuations. 



I. INTRODUCTION 

The large angle bump in the polarization-temperature 
cross correlation power spectrum [1] measured in WMAP 
data [2] indicates an optical depth to electron scatter- 
ing of 0.17 ± 0.04 [3]. To explain both this high opti- 
cal depth and the Lyman-a optical depth from Gunn- 
Peterson troughs [4] towards z ~ 6 quasars in the Sloan 
Digital Sky Survey [5] requires a complex reionization 
history [6-8]. While slight modifications to the large an- 
gular scale polarization power spectra exist with different 
reionization histories that integrate to the same electron 
scattering optical depth, one cannot use these changes 
to fully reconstruct the reionization history as a function 
of redshift due to large cosmic variance associated with 
anisotropy measurements at a few tens degree angular 
scales [9]. 

Under standard expectations for reionization, mainly 
due to UV light emitted by first luminous objects, the 
reionization process is expected to be both patchy and 
inhomogeneous [10]. This leads to fluctuations in the 
electron scattering optical depth and to a modulation 
of both temperature and polarization contributions such 
that new anisotropy fluctuations are generated at ar- 
cminute scales corresponding to inhomogeneities in the 
visibility function [11,12]. The increase in sensitivity and 
angular resolution of upcoming CMB anisotropy data 
suggests that such small-scale fluctuations can soon be 
used to understand the reionization history and associ- 
ated physics beyond measurements related to large-scale 
polarization alone. 

The secondary reionization-related anisotropies, unfor- 
tunately, are expected to be confused with other sec- 



ondary effects such as the Sunyaev-Zel'dovich (SZ; [13]) 
contribution at late times associated with galaxy clusters 
and first supernovae [14], and higher order effects such as 
gravitational lensing [15]. Though this confusion can be 
partly removed, such as through frequency cleaning in 
the case of SZ [16] or subtraction of lensing via higher 
order statistics [17], secondary anisotropies alone cannot 
be used to extract the detailed history of reionization 
beyond the integrated optical depth [12,18]. 

An interesting possibility involves the cross-correlation 
between small-scale CMB polarization maps and images 
of the high redshift universe. This cross-correlation can 
help distinguish broad aspects such as whether the uni- 
verse reionized once or twice [23]. Beyond the reioniza- 
tion history, mainly in the scattering visibility function 
as a function of redshift, it would be interesting to study 
additional details related to the reionization process such 
as the relative distribution of free-electrons and neutral 
Hydrogen. Since the brightness temperature fluctuations 
associated with the 21 cm spin-flip transition [19] trace 
the neutral Hydrogen distribution, the combination in- 
volving secondary CMB and 21 cm fluctuations [20] could 
then provide additional details related to reionization. 

In this paper, we consider such a combined study in 
the form of a cross-correlation between arcminute scale 
CMB anisotropy maps and images of the high redshift 
universe around the era of reionization from the 21 cm 
background. We make the assumption that the reion- 
ization process is inhomogeneous and patchy such that 
at certain epochs during partial ionization, contributions 
are generated to both the 21 cm background and CMB. 
In such an era, we also assume that the spin temperature 
of the neutral Hydrogen distribution is decoupled from 
CMB and is dominated by the kinetic temperature; In 
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this case, the signature in 21 cm emission is fluctuations 
related to the distribution of neutral Hydrogen, though, 
in the limit where the spin temperature is smaller than 
that of CMB (prior to reionization and appearance of 
first sources), 21 cm fluctuations will be tightly coupled 
to that of large angular scale CMB temperature; In this 
case, one naturally expects a perfect cross-correlation be- 
tween these two maps. We do not consider this possibility 
as the reionization process is expected to heat the IGM 
by z ~ 20. Also, opportunities for very low frequency ob- 
servations, where 21 cm signatures from rcdshifts prior 
to reionization are expected, are extremely limited. 

In the era of partial reionization, the existence of the 
proposed cross-correlation is due to the fact that inhomo- 
geneities that lead to fluctuations in both CMB (in terms 
of the electron distribution) and 21 cm background (via 
the neutral content) trace the same underlying density 
field. One expects the two to spatially correlate though 
this would be an anti-correlation as regions containing 
free-electrons will be mostly free of neutral material. This 
fact is captured by a spatial cross-power spectrum be- 
tween free-electron and neutral Hydrogen fluctuations. 
Since the 21 cm background allows one to probe the 
power spectrum of neutral Hydrogen alone, while CMB 
probes the power spectrum of reionized patches, in com- 
bination, the cross-power spectrum provides additional 
information on physics related to reionization. 

Unfortunately, while there is a strong (anti) cross- 
correlation between free-electrons and neutral Hydrogen 
fluctuations, we find the observable projected angular 
cross-spectrum to be insignificant due to a geometric can- 
cellation. This cancellation comes from the CMB side 
and involves the line of sight projection of the velocity 
field during reionization. The same cross-correlation can 
also be studied using a higher order correlation in the 
form of a bispectrum in Fourier space involving the large 
scale velocity field of ionized regions from the Doppler 
effect, arcminute scale CMB anisotropies, and the 21 
cm background. This measurement avoids the geometric 
cancellation associated with the line of sight projection, 
but due to the large cosmic variance associated with the 
velocity field traced by the Doppler effect, its measure- 
ment is limited to signal-to-noise ratios of order ten (with 
noise dominated maps) or, at most, a hundred. 

We discuss the measurement of proposed cross- 
correlations using CMB maps from upcoming missions, 
such as Planck surveyor and the South Pole Telescope 
(SPT), and maps of the 21cm background with, say, the 
Square Kilometer Array. Unfortunately, with signal-to- 
noise ratios around ten or below, it is unlikely that one 
can use the proposed cross-correlation bispectrum to eas- 
ily extract detailed information on the relative distribu- 
tion between neutral Hydrogen and electrons. 

The paper is organized as follows. In § II, we derive 
the existence of the cross-correlation both in terms of the 
cross power spectrum between CMB temperature and 
21 cm fluctuations and a higher order bispectrum that 
avoids a geometric cancellation associated with the cross 



power spectrum. In § III, wc discuss our results and 
suggest that though there is adequate signal-to-noise to 
perform a cross correlation study in terms of the bispec- 
trum. We briefly discuss how this measurement can be 
improved and what information related to reionization 
can be extracted from this measurement. We conclude 
with a summary in § IV. 



II. CALCULATION METHOD 

We first discuss fluctuations in the 21 cm background 
and then small-scale CMB anisotropies related to reion- 
ization. We will then consider the relation between these 
two quantities in terms of the angular power spectrum 
related to the cross-correlation. Beyond the power spec- 
trum, we also discuss a bispectrum associated with these 
two quantities and the large angular scale CMB fluctua- 
tions related to velocity variations during reionization. 



A. 21 cm Fluctuations 

In the low optical depth limit of the radiation transfer, 
the brightness temperature fluctuation associated with 
21 cm background can be written as [28] 
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where Tqmb = 2.73(1 + z) K is the CMB temperature at 
redshift z, Ts is the spin temperature of the IGM, and 
the optical depth in the 21 cm hyperfine transition, t Hi 
is 
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Here tiri is the neutral hydrogen density expressed in 
terms of fluctuations in the density as rim — XHn g {l + 
5 g ), where n g is the mean number density of cos- 
mic baryons and xr is the ionization fraction, vq = 
1420.4 MHz is the rest-frame hyperfine transition fre- 
quency, Aio is the spontaneous emission coefficient for 
the transition (2.85 x 10~ 15 s _1 ). We refer the reader to 
Rcf. [28] for further details. 

We expand fluctuations in the brightness temperature 
related to the 21 cm emission and write multipole mo- 
ments as 
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where xh is the mean neutral fraction. Making use of 
the Rayleigh expansion 
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we can now write the multipole moments as 
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where the window function is 
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In numerical calculations, we consider 21 cm measure- 
ments during partial rcionization (0 < xjj < 1) where 
the spin temperature is expected to be dominated by the 
kinetic temperature. In this case T s <~ T^ ^> Tcmb and 
the window function related to the 21 cm fluctuations is 
independent of the spin temperature of neutral Hydro- 
gen. 

In Eq. 6, r(z) is the conformal distance (or lookback 
time) from the observer at redshift z = 0, given by 
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where the expansion rate for adiabatic CDM cosmological 
models with a cosmological constant is 

H 2 = H 2 [O ro (l + zf + fl K (l + zf + fi A ] , (9) 

where H can be written as the inverse Hubble distance 
today Hq 1 = 2997.9/i _1 Mpc. We follow the conventions 
that in units of the critical density 3H 2 /8irG 7 the con- 
tribution of each component is denoted fit, i = c for the 
CDM, b for the baryons, A for the cosmological constant. 
We also define the auxiliary quantities TO = Q c + fib and 
flu = 1 — Y2i &i ! which represent the matter density and 
the contribution of spatial curvature to the expansion 
rate respectively. Although we maintain generality in all 
derivations, we illustrate our results with the currently 
favored ACDM cosmological model. The parameters for 
this model arc fl c = 0.30, fib = 0.05, fl\ = 0.65 and 
h = 0.7. 

Using the multipole moments, note that the angular 
power spectrum related to the 21 cm fluctuations is given 
by 
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such that 
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where P gg {k) is the power spectrum of neutral Hydrogen 
distribution. In Fig. 1(b), we show, as an example, the 
angular power spectrum of 21 cm fluctuations under the 
assumption that the neutral Hydrogen fluctuations trace 
the linear density field with a bias factor given by the halo 



based calculations [29] such that 
and the bias follows from 
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Here, M_ and M + are the lower and upper limits of 
masses, dN/dM is the mass function [21] and b(M, z) is 
the halo bias of Ref. [22] . For neutral Hydrogen fluctua- 
tions, we set M_ — > and M + the value corresponding 
to virial temperatures of 10 4 K. Note that we have taken 
a simple description of IGM bias here and the situation 
is likely to be more complicated. Since our main objec- 
tive is to see if there is an adequate signal-to-noise for a 
detection of the cross correlation between CMB and 21 
cm fluctuations, we ignore subtle details associated with 
neutral Hydrogen distribution. In Fig. 1(b), for refer- 
ence, we also show other contributions to the brightness 
temperature fluctuations at these low frequencies, such 
as CMB and free-free emission from ionized halos [30] . 

For cross-correlation purposes with secondary order ef- 
fects in CMB (discussed in the next Section), one needs 
to consider further fluctuations in the 21 cm background. 
In the patchy rcionization scenario, these are naturally 
present during the epoch when the universe transits from 
a neutral one to a fully reionized one. Making use of the 
inhomogencitics in the neutral gas distribution, we write 
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in terms of fluctuations in the neutral Hydrogen fraction. 
The brightness temperature fluctuations in the 21 cm 
background are 
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with the same window function as before (and again in- 
dependent of the spin temperature). In Eq. 13, the multi- 
plication between the density fluctuations and ionization 
fraction inhomogencitics in real space has been converted 
to a convolution between the two fields in Fourier space. 
Setting, ki = k— k' and k2 = k', the multipole moments 
can be simplified as 
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where we have defined a general integral over spherical 
harmonics written such that 

KC 2 u mi ^)= J dnY^{n)Y^{n)....Y^{h). (15) 



B. CMB fluctuations 

The bulk flow of electrons that scatter the CMB pho- 
tons, in the reionized epoch, lead to temperature fluctu- 
ations through the well known Dopplcr effect [24] 
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where v g is the baryon velocity. The visibility function, 
or the probability of scattering within dr of r, is 
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Here r(r) = drr is the optical depth out to r, X e (z) 
is the ionization fraction, as a function of redshift, and 
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is the optical depth to Thomson scattering to the Hubble 
distance today, assuming full hydrogen ionization with 
primordial helium fraction of Y p (= 0.24). 

Here, we make use of a description for X e (z) that is 
consistent with current data such as the WMAP optical 
depth [3], following Rcf. [7]. The description is based 
on the Press-Schechter [21] mass function and related to 
the reionization by UV light from the first-star forma- 
tion. Here, X e {z) varies from a value less than 10 _1 at 
a redshift of 30 to a value of unity when the universe 
is fully reionized, at a redshift of ~ 5. The reionization 
history is rather broad and the universe does not become 
fully reionized till a late time, though the reionization 
process began at a much higher epoch. The total optical 
depth related to this ionization history is ~ 0.15 consis- 
tent with WMAP analysis. The scattering is distributed 
widely given the long reionization process. This is the 
same model as model A of Rcf. [23]. Given the ioniza- 
tion fraction, we set the neutral fraction, related to the 
21 cm background, to be simply xh(z) = 1 — X e (z). 

In Fig. 1(a), wc show the linear Doppler effect, includ- 
ing contributions resulting due to double scattering effect 
described in [24] (see, Ref. [25] for details) following the 
reionization history. The power spectrum is such that it 
peaks around the horizon at the reionization projected on 
the sky today. As shown in Fig. 1(a), the Doppler effect 
cancels out significantly at scales smaller than the hori- 
zon at reionization surface since photons scatter against 



the crests and troughs of the perturbation resulting in a 
net cancellation of the anisotropy. 

The patchy reionization temperature fluctuations can 
be written as a product of the line of sight velocity, under 
linear theory, and fluctuations in the electron fraction, x e 
as 
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where we have relate the velocity field to the density field 
using the continuity equation such that 
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where the overdot represent the derivative with respect 
to radial distance r and the linear density field is scaled 
as a function of time using the growth function G(z), 
where S(k, r) = G{r)S(k, 0) [27] 
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Similar to the Ostriker-Vishniac (OV) effect [26], this 
second order effect related to reionization, called the 
patchy-reionization contribution, avoids the strong can- 
cellation associated with the linear Doppler effect. As 
discussed in Ref. [12], for high optical depth reioniza- 
tion histories, where the reionization process is expected 
to be highly inhomogencous, this effect dominates sec- 
ondary anisotropy fluctuations after the SZ contribution 
related to galaxy clusters; it is higher than both the OV 
and Kinetic SZ contributions. We show these contribu- 
tions in Fig 1(a). When calculating the patchy reion- 
ization power spectrum, we model the free-electron dis- 
tribution with a power spectrum that involves the same 
halo-based bias as in Ref. [12]; the bias is dominated by 
halos with temperature at the level of 10 4 K and above, 
where atomic cooling is expected and first objects form 
and subsequently reionize the universe. Since what en- 
ters in the power spectrum is the bias factor times the 
growth of density perturbations, in this case the product 
bx c (z)G(z) is in fact a constant as a function of redshift 
out to z of 20 or more given the rareness of massive halos 
at high rcdshifts [14]. 

We can expand the temperature perturbation due to 
reionization related scattering, T Rel , into spherical har- 
monics given by ai m — J dnr(n)Fj m *(n) such that 
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where we have symmetrizised by using ki and k2 to rep- 
resent k — k' and k' respectively. Using 
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we can further simplify and rewrite the multipole mo- 
ments as 
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C. Cross Power Spectrum 

We can now construct the angular power spectra by 
considering the clustering aspects between neutral Hy- 
drogen in the 21 cm background and fluctuations in the 
secondary reionization related effects. Under the assump- 
tion that the two fields are statistically isotropic, the 
cross-correlation is independent of m, and we can write 
the the cross-power spectrum between 21 cm fluctuations 
and the secondary CMB anisotropy fluctuations as 
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The cross-correlation can be written using 
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We can separate out the contributions such that the to- 
tal is made of correlations following (v g 5s)(8 Xe S XH ) and 
(v g S XH ){SsS Xc ) depending on whether we consider cumu- 
lants by combining ki with k^ or k 2 respectively. After 
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FIG. 1. (a) Power spectrum for the temperature 
anisotropics in the fiducial ACDM model with r = 0.15. 
In addition to primary CMB fluctuations, we also show the 
Doppler, SZ and Patchy Reionization contributions at late 
times during the reionization epoch and afterward. For ref- 
erence, we also show the instrumental noise power spectra 
of Planck and the South Pole Telescope; We will use these 
for signal-to-noise calculations, (b) The brightness temper- 
ature fluctuations at z/ bs = 140 MHz. Here, we show the 
21 cm signal in a 1 MHz bandwidth window around z — 10, 
the CMB, and expected foreground contamination related to 
free-free fluctuations in ionized halos at redshifts around 3 
and below. For reference, here we also show the instrumental 
noise contribution related to Square Kilometer Array (SKA) 
observations of the 21 cm fluctuations. 



some straightforward but tedious algebra, and noting the 
orthonormality of the Wigner-3j symbols 
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we can write 
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Here, we have only considered the contribution related to 
the (v g S}(x e XH) term; we ignore the contribution related 
to the {v g Xn){5x e ) term due to the mismatch in fluctua- 
tion scales between velocity fluctuations (that trace large 
angular linear scales) and the distribution of neutral Hy- 
drogen (which is expected to be patchy and concentrated 
to smaller halos); By ignoring this term, we have underes- 
timated the cross-power by, at most, a factor of 2 and this 
difference is unlikely to change our conclusions regarding 
the detectability. In simplifying the integrals involving 
spherical harmonics, we have made use of properties of 
Clebsh-Gordon coefficients, in particular, those involv- 
ing 1 = 1. The integral involves two distances and two 
Fourier modes and is summed over the Wigner-3j symbol 
to obtain the observable cross-power spectrum. 

Similar to the Limber approximation [31], in order to 
simplify the calculation, we use an equation involving 
completeness of spherical Bessel functions: 
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where the assumption is that F(k) is a slowly-varying 
function. Applying this to the integral over k 2 gives 
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where the mode-coupling integral is 
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D. Bispectrum 

As we will soon discuss, the cross power spectrum is 
not an ideal probe of the cross-correlation between fluc- 
tuations in the ionized and neutral regions around the 
era of reionization due to the geometric cancellation as- 
sociated with large-scale velocity projection. Instead, we 
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FIG. 2. The integrals J t (r), top lines (Eq. 31), and L t (r), 
bottom lines (Eq. 40), as a function of the multipole when 
r is given by the redshift indicated. The mode coupling 
integral related to the cross-power spectrum, Ji(r), oscil- 
late between positive and negative values and we label posi- 
tive/negative parts with (+)/(-). These functions captures 
the cross-correlation between velocity field related to the 
CMB secondary anisotropies and either the density field fluc- 
tuations, in the case of the cross-power spectrum, or with the 
velocity field traced by the large scale Doppler contribution, 
in the case of the bispectrum. 

consider the possibility for a measurement of the three- 
point correlation function or, in Fourier space, the bis- 
pectrum that is constructed in a way such that it avoids 
the cancellation related to the cross power spectrum. For 
reference, the three-point correlation function is related 
to multipole moments via 

B(n,m,i) = <T(n)T(m)T(i)) (32) 
ee ^(a limiffl|2m2ffii3m3 )^(n)^(m)^(l) , 

where the sum is over (Zi,mi), (Z 2 ,m 2 ), (h,ms). Statis- 
tical isotropy again allows us to express the correlation 
in terms an m-independent function, 



h h h 

mi 1Tl 2 iti 3 



Bi ± i 



(33) 



where the quantity -Bi 1 z 2 > 3 is described as the angular 
averaged bispectrum [25] and we have made use of the 
orthonormality relation of Wigncr-3j symbols (Eq. 27). 
As discussed in Ref. [25]. the angular bispectrum, Bi t i 2 i 31 
contains all the information available in the three-point 
correlation function and frequently used quantities such 
as the skewness and the collapsed three-point function 
can be written in terms of the bispectrum (such as a 
filtered sum of certain bispectrum configurations) [33]. 

We consider a bispectrum of the form Doppler- 
Rcionization-21 cm fluctuations. The bispectrum can be 
written as 
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1 



IT dop (fci,r) = -rr-gG. 



(36) 



In simplifying the integrals involving spherical harmon- 
ics, note that we have again made use of the properties of 
Clebsch- Gordon coefficients, in particular, those involv- 
ing 1 = 1. 

The bispectrum is given by 



E 



h h h \ /„21cm„dop ^Rei 
TTli m 2 m 3 



a hm 1 a l 2 m 2 a l 3 m 3 



(2Z X + 1)(2Z 2 + 1)(2Z 3 + 1) (h h h , , 

0' hM ' 



47T 



(37) 



Here, 



b h,h = ^2 J k i dk i J kidk 2 P88(ki)Px H x c ( k 2) 

x/«X(fci,fc 2 )/ i 2 1 lcm (fci)/ ; d 2 ° p (fc 2 ). (38) 

As with the cross-power spectrum calculation, we use 
Limber approximation to describe the coupling between 
21 cm and rcionization ionized and neutral fractions. The 
bispectrum term is then, 



II 



dr 



(gGG)W 21cm 



d\ 



dA ' 



; r 



Li 2 (r) 



(39) 



where the new mode coupling integral is now 

J' h (kr) f . j' h {kn) 



L h (r)= J k 2 dkP 5S (k) J -^l J dnigGXny- 



k 



(40) 



This can be directly compared to the mode coupling in- 
tegral associated with the cross-power spectrum. 




FIG. 3. The cumulative signal-to-noise ratio as a function 
of the multipole I. With thin lines, we show the case for the 
angular cross power spectrum between 21 cm fluctuations and 
small-scale CMB anisotropies while with thick lines, we con- 
sider the case involving the bispectrum between 21 cm fluc- 
tuations, small-scale CMB, and large-scale CMB (in this case 
the multipole shown in the x-axis is that of 1$ with the two 
other multipoles summed over). In both cases, the top lines 
are for a case involving all-sky observations with no instru- 
mental noise, but limited by the cosmic variance. The middle 
long-dashed lines assume instrumental noise level of the South 
Pole Telescope (SPT), shown in Fig. 1(a), combined with SKA 
noise for over 4000 sqr. degrees. The bottom lines show the 
case involving Planck data when combined with a map of the 
SKA 5 sqr. degree field of view. The signal-to-noise ratio for 
the cross power spectrum is smaller than that of the bispec- 
trum. 



III. RESULTS 

In Fig. 1, we show power spectra of CMB temperature 
anisotropies (top panel) and brightness temperature fluc- 
tuation spectra at 140 MHz corresponding to rcdshiftcd 
21 cm line emission from a redshift ~ 9 (bottom panel). 
When calculating 21 cm fluctuation spectrum, we assume 
a bandwidth for observations of 1 MHz centered around 
140 MHz; this roughly corresponds to a redshift inter- 
val of 0.06. For reference, in Fig 1(b), we also show the 
noise level related to a single month of observation with 
the planned Square Kilometer Array (SKA); We take the 
approach of Ref. [28] and assume a system temperature 
of T sys = 200 K. 

In Fig. 2, we illustrate the mode coupling integrals re- 
lated to both the cross power spectrum and bispectrum 
between 21 cm fluctuations and CMB anisotropies re- 
lated to reionization. In the case of the cross power spec- 
trum, the mode coupling term (Eq. 31) involves an inte- 
gral over the product of j[ and ji terms. The j{ term here 
is related to the line of sight projection of velocity fluctu- 
ations from the CMB side while the ji term is associated 
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with the line of sight projection of density fluctuations 
in the 21 cm map. The two projections, unfortunately, 
are not well correlated (analytically, this can be under- 
stood based on the phase shift between j[ and ji terms, 
while, physically, the density inhomogeneities peak at a 
smaller physical scale than the velocity field). The re- 
sulting effect is to produce a mode coupling integral that 
oscillate between positive and negative values. Since the 
final cross power spectrum is a weighted summation over 
this mode coupling term, this results in a net cancella- 
tion and a reduction of the cross correlation. Instead 
of the all-sky derivation used here, the same calculation 
can be considered under a flat-sky approximation. In this 
case, the correlation term involves a term that scales as 
/ d 3 k J dh([h-v\S(k)), which simplifies at the small angle 
limit to a zero when one integrates over the cos(#) angle 
between the line of sight and the velocity field. 

The relative cancellation between velocity and den- 
sity fluctuations can be avoided if one introduces an 
additional probe of the velocity field to the correla- 
tion measurement. We have suggested the bispectrum 
where, in addition to 21 cm fluctuations and small scale 
CMB anisotropics related to reionization, one makes use 
of large angular scale CMB anisotropies related to the 
Doppler effect. The mode coupling integral now involves 
two j[ terms that capture the velocity field from both the 
large scale Doppler effect and small scale anisotropies 
which are due to a modulation of the same velocity 
field. The mode coupling integral no longer oscillates and 
avoids the cancellation related to the cross power spec- 
trum. While this is an advantage, there is a disadvantage 
related to this bispectrum since it involves a map of large 
scale temperature fluctuations in addition to small scale 
anisotropies. In terms of the detection, the signal-to- 
noise ratio is limited by the cosmic variance related to 
the large scale Doppler effect which can be understood 
mainly as a confusion between these fluctuations and the 
ones generated at degree angular scales at the surface of 
last scattering such as the acoustic peak structure. 

In order to consider the extent to which the cross- 
correlation can be detected and studied to understand 
the reionization process, we estimate the signal-to-noise 
ratio for a detection of both the cross power spectrum 
and the bispectrum. In the case of the power spectrum, 
the signal-to-noise ratio is 



C, 



,Lot 



£ii _|_ ^moisc _|_ ^fforeg 



(43) 



;(2/+i)/ sky 



21cm-Rci^ 2 



-/force 



^21cm-Roi^ 2 _|_ £rCMB,tot£r21cm,tot 

(41) 



while the signal-to-noise ratio for the bispectrum is cal- 
culated as 

2 / n \2 



sky 2^ 



^21cm,tot£jCMB,tot£,CMB,tot ' 
WI2I3 h h h 



(42) 



Here, C\ represents all contributions to the power 
spectrum of the ith field, 



confusing foreground contribution, and / s k y is the frac- 
tion of sky covered by the data. In each of the cases 
involving CMB or 21 cm maps, we make use of the noise 
and foreground spectra shown in Fig. 1; In the case of 
CMB foregrounds, we make the assumption that SZ con- 
tribution can be perfectly removed based on multifrc- 
quency maps (such as with Planck). This is important 
as this leads to an improvement in the signal-to-noise ra- 
tio at arcminute angular scales with a reduction in an 
additional confusion. 

The signal-to-noise ratios are summarized in Fig. 3 for 
both the cross power spectrum and the bispectrum. Here, 
in addition to realistic scenarios, we also consider the 
maximal signal-to-noise ratio one could achieve by us- 
ing no instrumental noise all-sky maps; In the case of the 
cross power spectrum, the maximum signal-to-noise ratio 
is order a few while, for the bispectrum, this is at a level 
slightly below 100. In addition to the partial reduction 
of the geometric cancellation associated with the mode 
coupling integral, the signal-to-noise ratio with the bis- 
pectrum is higher as one has more modes from which to 
construct and measure the bispectrum when compared 
to that of the power spectrum. 

Using Planck and SPT maps, combined with SKA 21 
cm maps, we also show the level of the expected signal-to- 
noise ratios in Fig. 3. The bispectrum constructed with 
the 4000 sqr. degree map of SPT, combined with Planck 
large scale CMB map, could provide the first opportunity 
to detect the suggested cross-correlation. We note here 
that techniques now exist to construct the bispectrum re- 
liably from CMB and large scale structure maps and they 
have been successfully applied to understand the presence 
of non-Gaussianities in current data. Thus, we do not 
consider the measurement of the proposed bispectrum, 
involving CMB temperature (both at large and small an- 
gular scales) and the 21 cm background to be any more 
complicated than what is already achieved [32]. While 
the measurement is not complicated, unfortunately, since 
the cumulative signal-to-noise ratio is at the level of ten, 
it is unlikely that a detection of the bispectrum can be 
used to establish in detail the cross spatial power spec- 
trum of neutral and ionized regions during reionization. 

While the signal-to-noise is not high, we still recom- 
mend the proposed measurement in future data as we 
'may have underestimated the cross-correlation, and thus 
the expected signal-to-noise ratio, by a variety of effects. 
For example, if there is a significant contribution to CMB 
anisotropies from reionization, beyond the patchy reion- 
ization contribution we have suggested, such as in the SZ 
map, then this would be evident in the cross-correlation. 
The cross-correlation can also be used as a mechanism 
to separate out small scale CMB anisotropies related to, 
for example, lensing effect on CMB from that of reion- 
ization; this separation is generally hard to achieve since 
both effects have the same thermal spectrum and peak 
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at similar angular scales. 



IV. SUMMARY 

During the transition from a neutral to a fully reion- 
ized universe, scattering of cosmic microwave back- 
ground (CMB) photons via free-electrons lead to a new 
anisotropy contribution to the temperature distribution. 
If the reionization process is inhomogeneous and patchy, 
the era of reionization is also visible via brightness tem- 
perature fluctuations in the redshifted 21 cm line emis- 
sion of neutral Hydrogen. Since regions containing elec- 
trons and neutral Hydrogen are expected to trace the 
same underlying density field, the two are (anti) corre- 
lated and this is expected to be reflected in the anisotropy 
maps in terms of a cross-correlation between arcminutc- 
scale CMB temperature and the 21 cm background. 

In terms of the angular cross-power spectrum, unfor- 
tunately, this correlation is insignificant due to a geo- 
metric cancellation associated with second order CMB 
anisotropies. Thus, it is unlikely that the cross- 
correlation spectrum between small-scale CMB and 21 
cm fluctuations will be measurable even with maps in- 
volving a high signal-to-noise per pixel. The same cross- 
correlation between ionized and neutral regions, however, 
can be studied using a bispectrum involving large scale 
velocity field of ionized regions from the Doppler effect, 
arcminute scale CMB anisotropies involving reionization 
signal, and the 21 cm background. While the geomet- 
ric cancellation is partly avoided, the signal-to-noise ra- 
tio related to this bispectrum is reduced via large cosmic 
variance related to the velocity fluctuations traced by the 
Doppler effect. Unless velocity fluctuations can be inde- 
pendently established, it is unlikely that the correlation 
information related to the relative distribution of ionized 
electrons and regions containing neutral Hydrogen can 
be obtained with a combined study involving CMB and 
21 cm fluctuations. 
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